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PROTON RECOIL MEASUREMENTS OF THE PuBe SOURCE NEUTRON SPECTRA 
by R. L. Alexander , D. F. Shook, and C. H- Ford 
Lewis Research Center 

ABSTRACT 

The spectrum of a 5 Ci PuBe neutron source has been measured from 40 keV 
to 800 keV. The measurement was made using a spherical proportional counter 
containing 4 atmospheres of hydrogen* Pulses due to gamma-ray interactions 
in the chamber were distinguished by rise time analysis using a two-parameter 
pulse-height analyzer* The proton recoil spectrum was analyzed using a 
code written by Benjamin which accounts for contributions to the proton recoil 
spectrum from neutrons of energy greater than that to be analyzed* 

The proportional counter measurement is combined with a recent liquid 
scintillator measurement to give a complete spectrum. Good agreement exists 
between the two measurements in the region of overlap from 500 to 800 keV. 

The proportional counter data show that the flux is fairly smooth and flat 


in the region from 40 keV to 800 keV, 

The present measurements indicate 
source of 20.2 + 3 . 9 percent. 


INTRODUCTION 

For many applications of radioactive 
neutron sources it is necessary to have an 
accurate knowledge of the neutron spectrum 
as well as the source strength. The use of 
large water bath techniques has made it 
possible to determine source strengths to an 
accuracy of 3 percent or better but the 
neutron spectrum is not nearly so well known. 
Considerable past work with emulsions and 
organic scintillators has resulted in fairly , 
well established spectra at energies greater 
than about 1 MeV but below this energy the 
spectra of many radioactive sources have not 
been measured. 

In the present work a .5 Curie (Ci) PuBe 
source spectrum is measured from 800 keV 
down to 40 keV using a spherical proportional 
counter (l). This measurement is combined 
with a recent liquid scintillator measure- 
ment (2) to give a complete spectrum. The 
primary difficulties in making a proportional 
counter measurement of the low energy part 
of the Be(oC,n) spectrum are caused by the 
large number of high energy neutrons emitted 
by these sources, and by the gemma -ray ■ 
background* 

EXPERIMENT 

A spherical proportional counter de- 
signed by Benjamin et.&l*, (1) ms used in the 
measurements* The counter ms filled with 
4 atmospheres of hydrogen and a trace of 


a neutron yield below 1 MeV for this 


helium -3 an d had an inside diameter of 
3*94 cm* The counter wall was .051 cm thick 
stainless steel. A lead shield about . 64 cm 
thick, and an outer cadmium cover about 
.075 cm thick surrounded the counter. 

The measurements were made with the 
counter mounted on a light frame 1 . 8 meters 
above the floor of a room 4 meters high, 
and 6 meters square. The PuBe source was 
supported with its axis vertical at 15 cm 
from the counter by a . 62 cm diameter 
threaded iron rod* 

A block diagram of the electronics is 
shown in figure 1. To separate protons from 
electron recoils, a two-parameter analyzer 
was driven by an identification signal as 
well as by a linear signal* To preserve the 
linear response for large pulses in hydrogen, 
which may take up to 6 / 4 . sec to rise, the 
preamplifier decay time was made 130 fx sec, 
and a single 7 ju sec RC integration and 
differentiation time constant ms chosen 
in the linear amplifier. To obtain an 
identification pulse separating fast-rising 
proton recoils from electrons, the pre- 
amplifier output was R-C clipped with a short 
time constant in the second linear amplifier. 
This pulse was used to drive a fast stretcher 
and a slow stretcher of about 3 /a sec 
duration. The slow stretcher pulse was 
adjusted to reach its peak at its end* 
Finally, the identification pulse was delayed 
4* 75 sec* These precautions ensured that 
the two- paramet e r analyzer would always 
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convert the full voltage of the linear 
signal «, 

Careful measurements were made of the 
analyzer linearity and pulse height zero. 
Amplifier gains were measured and absolute 
gas. gains were obtained using the 764 keV 
helium-3 peak® A plot of In A/V versus In V 
gave a straight line over a range of gas gain 
from A=8 to A=250. 

Measurements were made with gas gains of 
from 45 to 180. The resolution of the 
helium-3 peak was about 10 percent full width 
at half maximum (FWHM) * 

TEE 5 CURIE PuBe SOURCE 

This source is a right cylinder, outer . 
diameter 3,33 cm, and length 6.91 cm* It 
contains 79,88 gms of plutonium and 39,30 gms 
of beryllium intimately mixed and melted in a 
tantalum cup and doubly encapsulated in 
stainless steel. 

The source emission rate was 9,50 + ,29 
neutrons per second in September of 1965, 

The isotopic fraction of plutonium-24l was . 

. 74 percent (3) . This corresponds to an 
initial rate of increase of source strength 
due to the formation of americium- 241 of 2,1 
percent per year. 

The source strength in September 1970 
is estimated to be 13. 6 percent greater than 
in 1962. The relative emission rate of this 4 
source as a function of direction was measured 
with a long counter, the center of the front 
fact of which was at 100 cm from the center 
of the source. The emission rate in the 
direction perpendicular to the source axis 
was measured to be I.O67 + *020 times as 
great as the emission rate averaged over the 
total solid angle. The calculated neutron 
flux at one meter from the source perpen- 
dicular to the axis is then 

A , 9: . 50xl. 061X1. iy»d j. _ 9l i 
t i'rr (lOO)^ 

MIA ANALYSIS 

The maximum pulse size that can be pro- 
duced by a gamma-ray in this chamber cor- 
responds to that of about a 100 keV proton, 
so that there is no gamma- ray interference 
above this energy. Consequently, the proton 
recoil data analysis is straightforward. 

The Spec 4 code of Benjamin (4) ms used. 

This code requires a neutron spectrum shape 


above the highest energy to be analyzed. 

The required shape ms obtained from the 
data of reference 2. The code generates a 
proton recoil response shape for neutrons 
above the energy to be analyzed, using 
analytic response functions. It then 
normalizes this shape to the data, and strips 
the differentiated remainder. 

Below about 100 keV proton energy, 
gamma events occur in the chamber. A two- 
dimensional plot of counts from a cobalt-60 
source Is shown in figure 2. The pulses 
from a 40 keV electron travel a total 
distance equivalent to the diameter of this 
chamber, and may require* as long as 6^ sec 
to rise. Consequently, at a given energy 
the clipping circuit places most of the 
pulses near the energy axis. 

Figure 3 shows data from the PuBe 
source. A 100 keV proton in the counter 
travels only about a millimeter, so that 
all of the 100 keV proton pulses should 
rise within the clipping time of the identi- 
fication circuit. These pulses appear as a 
peak in the identification channel. Since 
the identification pulse size for these 
pulses is proportional to the energy, this 
peak occurs in a diagonal line on the plot (5). 

Unfortunately, not all the pulses on 
this plot can be clearly ascribed to the 
proton-recoil peak, or to gamma- rays. 

Above 100 keV, where the gammas ray induced 
pulses are gone, there is a flat continuum 
of pulses which have smaller identification 
pulses than those in the peak# It is be- 
lieved that these pulses are due to higher 
energy protons that lose only a fraction of 
their energy in the gas before colliding 
with the wall. At 2.5 MeV for example, a 
proton will lose 100 keV while travelling 
nearly 2 cm, in the gas. 

■Figure 4 shows the spectrum of pulses 
from the source surrounded by a spherical 
shell of tungsten, which sharply reduces 
the fraction of high-energy neutrons in the 
spectrum. The calculated fraction of the 
flux above 2 MeY in the tungsten leakage 
spectrum Is only 17 percent compared with 
73 percent for the bare source. The number 
of counts not on the diagonal have been much 
reduced. This non-peak fraction is 35 per- 
cent at 100 keV for the bare source, and 
only 4 percent for the tungsten leakage 
spectrum* 
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It is possible to subtract an arbitrary 
fraction of a gamma- ray spectrum from a 
mixed spectrum so as to obtain a smooth con- 
tinuation below 100 keV of the distribution 
above kO keV; so that the subtraction is 
from the non-peak counts* 

Figure 5 shows the total counts 
measured for the 5 Ci PuBe source. Two 
slightly different g&msm fractions are sub- 
tracted below 100 keV. The non-peak con- 
tribution to the spectrum for these sub- 
tractions is shown also* A modified version 
of Spec 4 was used to calculate the fraction 
of recoils from neutrons of energy greater 
than 1.5 MeV, This is also shown in figure 
5. Figure 6 shows similar results for the 
tungsten leakage spectrum; where the non- 
peak fraction is much lower. 

In both fluxes^ the fraction of non- 
peak, counts is not a constant with energy. 
For a fixed clipping tim®; there will be 
some energy of proton- recoil below which all 
the recoils will be in the group along the 
diagonal. The rise time for these very 
small pulses will be within the clipping 
time. Since the absolute number of recoils 
from high energy events increases only about 
10^ percent from 100 keV to 10 keV; the 
number of counts occurring off the diagonal 
should^ at most; rise no more than this* 

The uncertainty in the gamm-ray sub- 
traction above 40 k@V is of the order of the 
difference between the two subtractions 
shown in figures 5 and 6, The approximate 
upper limit of the counts in the non-peak 
fraction is the upper curve shown. If one 
subtracts more than enough to give the lower 
curve; portions of the identification 
channel region become negative, 

RESULTS AMD DISCUSSIOM 

The solid line of figure 7 shows the 
neutron flux spectrum reduced to 1 meter 
distance for the 5 0± PuBe source as pre- 
sently measured by the 4 atmosphere hydro- 
gen proportional counter. The curve drawn 
through the data is the weighted averages 
of two runs each above and below about 100 
keV* The hatched area represents the esti^ 
mated error* The analysis is- not inconsis- 
tent with a nearly featureless spectrum. 

But the analyses of runs taken on different 
days show the source spectra features re- 
produced here. Below ?0 keV; the large un- 
certainty in the spectrum is caused by the 
uncertainty in the fraction subtracted 

from the non-peak data* Because the non* 


peak proton recoil contribution is expected 
to decrease with decreasing neutron energy $ 
the upper limit of the spectrum shown is 
believed to be near the upper limit for the 
spectrum in this region. The lower limit 
is not so clearly defined. 

The’ integral of the flux from 0 to 800 
keV with a flat 30 n/cjrpsec MeV below 40 
keV is 15.6 + 2. 7 n/cnrsec* The total flux 
at one meter”" is calculated to be 91* 6+3, 7 
n/cmsec, The data of reference 2; there- 
for©; have been normalized to include 
76,0 n/cm sec between 800 keV and 11 MeV. 

The fraction of the neutron spectrum below 
1 MeV is then 20,2+3*9 percent. 

The spectrum of reference 2 is very 
similar to that of Anderson and Bond (6) 
who measured the spectrum of a similar PuBe 
source with nuclear emulsions. Their 
measured spectrum extrapolated smoothly to 
zero neutron energy indicated about 18 per- 
cent below 1 MeV, St. Remain; et, al, ; (7) 
using Bonner spheres to estimate neutron 
yields found about 17 percent of the spec- 
trum of a 1 Ci PuBe source below X MeV. 

They inferred a peak in the neutron spectrum 
at 0*3 MeV due to inelastic scattering of 
alpha particles from beryllium. The excited 
beryllium nucleus was postulated to decay 
by neutron emission. Shook; et.al. ; (8) 
measured the age of neutrons in various 
media and using a similar 5 Cl PuBe source; 
estimated that 21 percent of the neutrons 
were emitted below 1 MeV, They used the 
Anderson and Bond spectrum above 1 MeV, 

In conclusion; the neutron spectrum of 
a 5 Ci PuBe source measured from 40 keV 
to 800 keV is generally flat and smooth. 

The measurement indicates there are about 
20,2 & 3® 9 percent of the neutrons emitted 
below” 1 MeV; which is consistent with other 
measurements in magnitude but not in shape. 
The method of estimation of the contribution 
of gamma- my s to the recoil spectrum below 
neutron energies of 100 keV is the principal 
source of uncertainty in the present proton- 
recoil measurements which have large 
fractions of the proton recoils caused by 
high-energy neutrons. 
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Figure 1. - Block diagram of instrumentation. 


IDENTIFICATION PULSE HEIGHT 


t 2 3*5 

1 1 3 1C 7 2 

2 185217 

3 18117* 

* 859*2 1 

5 *87*3 13 

6 336C6 1*9 5 1 

7 242*8 *95 *8 

8 18533 795 196 1* 1 


7 8 9 10 11 12 13 1* 15 16 17 18 19 20 21 22 23 2* 25 2'6 27 28 29 33 31 3 


9 13715 889 331 63 6 

10 10e2* 965 *25 125 J7 5 

11 8*25 993**0 250 60 18 3 

12 6*77 10*0 538 256 95 32 2 1 

13 S155 102* 556 250 H6 50 18 2 

14 *C53 1C07 585 33C 139 60 23 8 

15 3*7* 987 605 3*8 200 75 30 12 

16 2679 929 603 330 207 98 *5 1* 


40 


1 7 

18 

19 

20 
21 
22 
23 
2* 


2211 912 603 383 207 108 ** 20 9 3 1 
175C 867*00 3** 250 117 59 27 5 2 
133C 801 609 37C 238 11* 63 29 13 3 1 

1C9* 7*1 537 38* 218 130 64 35 17 6 5 1 

758 66* 526 329 209 128 82 2* 21 11 l 2 

577 563 *72 362 23* 138 77 32 22 tl 3 1 1 

* C* 511 **8 3*0 192 1*1 76 ** 25 5 6 l 

3C9 466 *30 269 209 117 78 50 31 10 5 1 


26 23* 33* 371 326 198 12* 7* 38 25 ll 5 2 l 

26 176 293 333 27* 199 110 7* 37 16 11 7 2 1 

27 91 222 236 22* 191 127 65.41 21 ll 8 * 1 

28 75 190 228 197 160 100 55 30 116 U 7 3 

29 55 152 196 169 132 10658 2* 18 9*11 

30 21 101 1*5 1*9 132 78 65 29 25 16 3 II 

33. i« 75 1*3 152 ll* 7* 51 35 17 U 3 * 3 

32 1C 66 115 lie 110 75 ** 27 19 75 5 l 


33 

3* 

35 

36 

37 

38 

39 
*0 


8 *0 87 96 90 66 *9 31 19 7 7 6 3 l 1 

1 20 63 68 80 60 *t 21 9 9 5 3 2 1 

l 22 *5 6C 65 5* 39 22 10 9 * l 1 1 

12 36 57 *9 *2 3* 2* 7 8 7 1 2 1 

1 10 2* *0 29 38 17 16 15 8 * 3 I 

l 2 1* 27 30 26 27 23 9 10 1 2 l 

4 12 l* 21 25 23 7 ll 7 5 2 2 2 2 

5 17 26 23 20 5 8 3 1 1 1 


*1 

*2 

*3 

*5 

*6 

*7 

*8 


15 9 

8 10 
6 12 
6 9 

5 6 

3 3 

2 

3 3 


8 IS 5 5 1 1 

8 10 * 2 1 1 

3 8*11 

8 6 16 2 

6 7 3 2 1 

* 2 3 2 1 

3 2 1 2 2 1 

2 2 3 3 1 
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*9 

50 

51 

52 

53 
5* 

55 

56 


* 2 
3 2 

1 3 
1 2 
3 l 

2 


2 3 

1 1 
1 2 2 
2 2 1 
2 

1 1 


124 


57 

58 

59 

60 
61 
62 
63 
6* 


12 III 
l 1 

1 111 
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Figure 2. - Proportional counter data, cobalt-60 source. 
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Figure 5. - Comparison of measured with calculated count rates for 5-Ci Pu-Be source. 
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